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By Yoichi Motohashi*) 



Abstract: The aim of the present article is to render the spectral theory of mean values of automorphic 
L- functions - in a unified fashion. This is an outcome of the investigation commenced with the parts XII and 
XIV, where a framework was laid on the basis of the theory of automorphic representations and a general 
approach to the mean values was envisaged. We restrict ourselves to the situation offered by the full modular 
group, solely for the sake of simplicity. Details and extensions are to be published elsewhere. 
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1. To begin with, we stress that all notations 
and conventions will stay effective once introduced. 
Let G = PSL(2, K) and T = PSL(2, Z). Write 



1 x 
1 



k[(9] 



cos 9 sin t 
- sin 9 cos ( 



and N = {n[x} : x € R}, A = {a[y] : y > 0}, K = 
{k[0] : 9 € M/ttZ} , so that G = NAK is the Iwasawa 
decomposition of the Lie group G. We read it as 
G9g = nak = n[a;]a[?/]k[0]. The coordinate (x,y,9) 
will retain this definition. The center of the universal 
enveloping algebra of G is the polynomial ring on the 
Casimir element Q — y 2 (d 2 + d 2 ) —yd x dg. The Haar 
measures on the groups N, A, K, G are normalized, 
respectively, by dn = dx, da = dy/y, dk = dO/n, 
dg = dnda.dk/ 'y, with Lebesgue measures dx, dy, d9. 

The space L 2 (r\G) is composed of all left J 7 - 
automorphic functions on G, vectors for short, which 
are square integrable over r\G against dg. Elements 
of G act unitarily on vectors from the right. We have 
the orthogonal decomposition 

L 2 (r\G) = C-10°i 2 (r\G)0 e i 2 (AG) 

into invariant subspaces. Here °L 2 is the cuspidal 
subspace spanned by vectors whose Fourier expan- 
sions with respect to the left action of N have van- 
ishing constant terms. The subspace e L 2 is spanned 
by integrals of Eisenstein series. Invariant subspaces 
of L 2 (r\G) and _T-automorphic representations of 
G are interchangeable concepts, and we shall refer 
to them in a mixed way. 



The cuspidal subspace decomposes into irredu- 
cible subspaces 



°L 2 (r\G) = ®v. 



The operator becomes a constant multiplication in 
each V so that Q\v°° — \y v — 4) T, where V°° is the 
set of all infinitely differentiable vectors in V. Un- 
der our present supposition, V belongs to either the 
unitary principal series or the discrete series; accord- 
ingly, we have Vy € iM. or vy = £ — | (1 < I € Z). 

The right action of K induces the decomposition 
of V into AT-irreducible subspaces 



V 



© V p . 

p— — oo 



dimVp < 1. 



If it is not trivial, V p is spanned by a r'-automorphic 
function on which the right translation by k[9] be- 
comes the multiplication by the factor exp(2ip(9). It 
is called a -T-automorphic form of spectral parameter 
vy and weight 2p. 

Let V be in the unitary principal series. Then 
dim V p = 1 for all p € Z and there exists a complete 
orthonormal system {ip p £ V p : p € Z} of V such 
that 



OO / N 

E ^B^ ssn(n V P (a[|n|]g;^), 



1— — 00 



where 4> p (g;v) = y v+ ^ exp(2ip9), and 

p OO 

A s <p p (g; v) = / exp(—2TTiSx)4>p(vfa[x]g;p)dx, 
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with w = k[^7r] the Weyl element. Our normaliza- 
tion is such that the coefficients gv(n) do not de- 
pend on the weight. Also we may impose the Hecke 
invariance; in particular Qv{— n) = evQv( n ) with 
e v = ±1. 

We have skipped the discrete series. In the se- 
quel as well, we shall argue as if there were no dis- 
crete series representations. This should not cause 
any confusion, as the discussions below, especially 
those pertaining estimations, extend readily to the 
discrete series. Nevertheless, it should be remarked 
that the definition of </? p (g, V) given above has to be 
modified for those V in the discrete series according 
to the normalization [2, (2.21) and (2.25)]. 

Our discussion depends much on the uniform 
bounds for A<fr p (a[y], v), A = A + , such as [2, (4.3) 
and (4.5)]. In order to make our argument applicable 
to any cuspidal representations, we derive from the 
latter a bound that is somewhat weaker than the 
former but is still sufficient for our purpose; in fact 
the proof of [2, (4.5)] works for all cases. We thus 
put 

T p (a,v)= f y'-^MplM,") — , Res>l; 

jo y 

we have shifted the argument in [2, (4.10)] by — ^. 
Let us assume that v e iR. We divide the integral 
at y = \p\ + \v\ + 1. To the part with smaller y 
we apply the fact that yt ssn ^0 p (a[|u|]; v) is a unit 
vector in L 2 (R X , d x /n), d x u = du/\u\ (see e.g., [2, 
Lemma 4]). Hence this part is < (\p\ + \v\ + l) Kcs ~ 1 . 
On the other hand, by [2, (4.5)] the remaining part is 
< (H + IH + l) Ros ~i We then invoke the identity 



T p (s, V) = 47T • 



■KT p {s + 2,v)-pT p {s + l,v) 



1)2 



which can be proved by applying the operator D„ = 
(d/dy) 2 - (2tt) 2 - (y 2 - \)y~ 2 either to y s or to 
A(f> p (a[y],h') in the last integral, on noting that the 
latter is a constant multiple of the Whittaker func- 
tion Wp.yiAny) (see [2, (2.16)]). Then, by Mellin's in- 
version, we conclude that A<p p (a,[y], v) <c y^~ £ (\p\ + 
\v\ + l) 2+e for any small e > 0. 

The above works with v = I — \, the discrete 
series, as well. We remark that 1> v is related to the 
operator dg via the Kirillov map. Hereafter we shall 
term [2, (4.5)] and the bound thus obtained the basic 
bounds. 



2. We now define the automorphic L-function 
associated with an irreducible T^-automorphic repre- 
sentation V by 



L v(s) = ev{n)i 



This converges absolutely for Re s > 1 and continues 
to an entire function which is of a polynomial order 
both in s and in vy if Re s is bounded. 

We fix an A among V's, and consider the mean 
square 



M(A,g) 



/OC 
\L(± + it)\ 2 g(t)dt, L = L A , 
-OC 



where the weight function g is assumed, for the sake 
of simplicity but without much loss of generality, 
to be even, entire, real on R, and of fast decay in 
any fixed horizontal strip. Our aim is to establish a 
full spectral decomposition of M(A,g). Our method 
is applicable to any representation, though we shall 
deal with only those A in the unitary principal series. 
With a minor modification, the fourth moment of the 
Riemann zeta-function could be discussed equally. 
We start as usual with the integral 



I(u,v;g) 



L(u + it)L(v + it)g(t)dt. 



This is an entire function of u, v; in particular 

M{A,g) = l(\,\;g). 
In the region of absolute convergence we have 
R(u + v) 



I(u,v;g) 



C(2(u + t>)) 



£(0) + J(u,v;g) + J(v,u;g), 



where R is the Rankin zeta-function attached to A, 
g the Fourier transform of g, and 



J(u,v;g) = ^^ {2n + f)u+v 

f=l n=l y J 1 

2a 



( y/n(n + /) 
^ 2n + f ) 

where g = qa and 

g(x; u, v) = 2 u+v+2a 



g(f/(2n + f);u,v), 



.9(log((I + x)/(I-x))) 
(l-x) u+a (l + x) v+a 
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with < x < 1. Here a is a sufficiently large posi- 
tive integer, which is implicit throughout the sequel. 
Note that we have slightly modified our procedure 
given in [6] . 

Let g* be the Mellin transform of g. It is im- 
mediate to see that g*(s; u, v) is of rapid decay with 
respect to s, provided Res and u, v are bounded; 
moreover, g*(s;u,v)/T(s) is entire over C 3 . Thus, 
by Mellin's inversion, 



J(u,v;g) 



-fit 



f s D f (u + v - s) > g*(s;u,v)ds, 



where (77) is the line Res = 77 > 0, and 



g{n)g(n + /) / y/n(n + J) 



^ (2n + /)« ^ 2n + f J 

Checking the convergence, we see that the condition 
Re (u + v) > max{2, 1 + 77} is required here. 

3. Now, we need to have a full spectral decom- 
position of Df(s) which is to yield a continuation, 
with a polynomial growth, to the domain Re s > vo 
with a vo < 1 so that J{u, v; g) can also be continued 
to a neighbourhood of the point (5,5)- Recently, V. 
Blomer and G. Harcos [1] succeeded in establishing 
such an assertion on Df(s). They started with an old 
idea of ours, published in [6], to employ the Kirillov 
model to pick up a favourable automorphic form in 
dealing with the problem of the spectral decomposi- 
tion of D/(s), and proceeded one step further by a 
use of the theory of the Sobolev norms with which 
they could derive the crucial polynomial growth that 
we had left open in [6]. 

In what follows, we shall take an alternative way 
by modifying [6] with certain simple devices from [2] , 
[5], and by reworking our relevant unpublished notes. 
The present article can be read independently of [1] . 

Thus, let A, a be as above, and r a parame- 
ter in the right half plane; all implicit constants in 
the sequel may depend on A, a and Rer at most. 
We apply the inverse Kirillov map to the function 
w(y,r) which is equal to y a+ 2 exp(— ry) for y > 0, 
and vanishes for y < 0. By [2, Lemma 4] or [6, 
Lemma 1], there exists an automorphic form $(g, r) 
in the subspace A such that 



$(n[ir]a[y],T) 



g(n) 



1 V™ 



w(ny, t) exp(2ninx). 



More precisely, 

00 

$ (g> T ) = E a p( r Vp(g'^)> 
p— — 00 

where 

a p( T ) = - I w a {y, T)A<j> p {a[y];v A ) —. 

* Jo y 

The function $(g, t) is regular for Rer > 0. In fact, 
we have 

a p (r)«(|r| + l) 2 "(H + l)-- 

and for any g it holds that (p p (g, A) <C + l) 2 , by 
[2, (4.5)] applied to A<p p (a[y]; va)- 

To prove the bound for a p (r), we use the oper- 
ator D„ again: we have, for p ^ 0, 



a p (r) 



1 



4tt 2 p J 



w a {y, T)D VA A(j)p{s\y\\ v A ) dy, 



Integrate in parts, and repeat the procedure as long 
as the resulting integrand vanishes at y = 0, that is, 
a times at least, and we get the bound. The case 
p = can be skipped obviously. 



4. Next, we put 5 ( (g, r) = $(g, r)$(g, r); this is 
regular for Rer > 0. We have, for any integer / > 0, 



Jo 



ty(n[x\&[y], t) exp(— 2irifx)dx 



= y 2a+1 J2e(n)Q(n + f)(n(n + f)r 

n=l 

■ exp(-(2n + /)ry). 

The Parseval formula implies that the left side is 
equal to 



J 

J(0 



V7 

1(0) V7C(1 + 2v) 4tti 

where V runs over all irreducible cuspidal represen- 
tations, and 

00 

A{y,r;V)= E T )' <Pr(-> V))A<j>My\; v v ), 

p— — OO 
OO 

A(y,T-,v)= E (y(;T),E p (;v))A<t> p (4y];v)- 
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Here ( , ) is the natural inner product on L 2 (T\G), with the Maass operator E = e™(2iyd x + 2yd y — 

ids), we see that Qipk^pi is a linear combination of 
(fk+jfi+j, j = — 1, 0, 1, the coefficients of which are 
polynomials of the second degree on k, I; note that 
va is now regarded as a constant. Thus 



and E p is the Eisenstein series of weight 2p (see [2, 
(3.19)]). 

The convergence of the spectral expansion of W 
is in fact absolute and fast, provided a is sufficiently 
large. To confirm this, we apply the operator tt + idg 
repeatedly as is done at [2, (5.9)]; and we get, via the 
above bound for a p (r), 



(*(.,r),M-.^)><(M + i) 4a (l«v| + bl)-' Q 

Then we appeal to the basic bounds, getting 



A(y,r;V) 

« (|r| + 1) 4Q (M + l)-i a yi- £ cxp(-^y). 

The contribution of the continuous spectrum or ra- 
ther the function A(y, r; v) is to be discussed later. 
This yields our assertion. Note that A(y,r; •) is reg- 
ular for Re t > 0. 

In order to pick up a particular Fourier coef- 
ficient, we have used the projection procedure as 
Blomer and Harcos did in [1] to avoid an inner prod- 
uct argument that was applied in [6]. It is worth 
pointing it out that in [2, Sections 3-5] a procedure 
of the same kind was developed, employing the Kir- 
illov model as a main implement to evaluate projec- 
tions of a certain Poincare series explicitly in terms 
of its seed function. Therefore, we surmise that the 
projections A(y, r; •) could also be handled more pre- 
cisely with a modification of the argument of [2] . To 
this issue we shall return elsewhere. 

5. We are now to make the last estimation pro- 
cedure explicit. Thus, we note that 



1 



(w 2 - \ - iC2 P yY 



<<(K| + b|)- 2 i(Q + ^)«*(-,T)|| 

for any integer q > 0. By definition 

oo oo 

fe=i 1=1 

with ip k (g) = ip k (g,A), tpi(g) = ipi{g,A). Since 



3 = -Q 



where dj(k, I; q) are polynomials of degree 2q on k, I; 
and 



k=l 1=1 

where 

<? 



3=-q 



«(|fc| + |*| + l) 2 « 



(M + i) 4 



((|fc| + l)(|Z| + !))«' 



for each fixed q. We then note as before that (p k (g) -C 
(|fe| + l) 2 ; and thus, with q = [|a], say, we get the 
uniform bound 

(Q + ld 2 )^(g,T)^(\T\+l) 4a , 

which gives the first inequality in the last section. 

6. As to the continuous spectrum, we invoke the 
functional equation for E p ([2, (3.33)]), and have, for 
Re v = 0, 

= lp (v) f {si + id 2 y*(g, t)e_ p ( s , v)d g , 

Jr\G 



with 



J p (v) = TV 



. 2l/ C(i + 2i/)r(| + i/+p) 

^\-2v)T{\-v+p)' 



Assuming Rev > \, we unfold the last integral, 
and see that it equals R(y + \)Y p {v, r; q)/((2u + 1), 
where R is as above, and 

oo 

Y p (v,T;q)= Y b i+p,i(T;q) 

l — — oo 

•V / A 5 (j)i +p {a.[y],v A )A s <i)i{a[y] 1 v A )y lJ -idy. 
8=± Jo 
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Again by the basic bounds, we see that Y p (y, r; q) is 
regular and <C (|r| + l) 4 " for Rct > and Rc^ > 
— 2- Hence, in the same domain, 

A(y,r;u) = ^±^ A(y,r ] u), 



with 



Y P {v, r; q) 



T{\ + v+p) 



T{\-u+p) 



A4> p {&[y\;v). 



One may conclude, via [2, (4.3) and (4.5)], that 

« (\t\ + i) ia (W\ + i)-i a y i-\*»"\-' exp(-Vy) 

for Rer > and Rez^ > — |. 

7. We now set r = s, and observe that 



Df(s) = 



„s+2a roc 



I' 

Jo 



T(s + 2a) 



j'OO 

/ V s - 2 

Jo 



^>(n[x]a[y],s) exp(2irifx)dxdy 



J2fi-* gv (f)Z(s,V) 



i (0) C(TTW^)^ + ^ (s '^' 



where 



^V) = f ^l° y^A(y,s-,V) d y, 



s+2a roo 



T(s + 2a) Jo 



y s z K{y,s;v)dy. 



The bound for A(y, r; V") implies that S(s, V) is reg- 
ular and < |s| 4a +i(|j/y| + l)-3« for Res > ±. Sim- 
ilarly, E(s,v) is regular and <C |s| 4a+ ^(|^| + 1)~5 Q 
for Res > \ + |Rct/| and Re;/ > — |. 

Therefore we have proved that Df(s) is indeed 
regular and of polynomial growth for Res > | ([1, 
Theorem 2]). 

It is to be observed that in order to offset the 
exponential growth of the factor \/Y{s + 2a) we have 



adopted an old idea of Ju.V. Linnik which he intro- 
duced in his investigation on approximate functional 
equations for Dirichlet L-functions. The introduc- 
tion of the parameter r was made with this effect in 
mind. 

8. We return to the function J(u,v;g); thus we 
impose Re (u + v) > max{2, 1 + 77} initially. In view 
of the fast decay of g*(s;u, v), the last assertion on 
Df(s) yields immediately that 

J{u,v;g) = ^L v (u + v - \) <d(u,v;V]g) 
v 

1 f ({u + v-± + v)Q(u + v-\-v) 



Atti 



(0) 



<(1 + 2i/)C(l - 2u) 
R{y + 5) @{u, v; v; g)dv, 



where 



Q(u,v;V;g) = ^- f S(« + v - £, V)g*(^, u, v)d£, 
2 " J( V ) 



and the factor E(u + v — £, v) appears instead in the 
continuous spectrum. The 9 is of fast decay either 
in Vy or in u, as S is. 

We fix a sufficiently small e > 0; we may move 
the last contour to (e), provided Re (u + v) > |, say. 
Hence, the last expansion of J(u, v; g) holds under 
Re(w + v) > |. This lower bound is required to 
get the factors Ly (u + v — |) and ((u + v — ^ ± v) . 
However, the former is entire and of a polynomial 
order in vy if u, v are bounded. Thus the cuspidal 
part of J(u, v; g) is regular in a neighbourhood of the 
point (i, |) at which it takes the value 



J2Ly(h)®A(V;g), 



with Q A {V;g) = Q{\,\;V;g). 

As we are about to deal with the continuous 
spectrum, we should remark that Q(u, v; v; g) re- 
mains regular in the three complex variables and of 
fast decay in v, throughout the procedure below, be- 
cause of the property of S(s, v) mentioned above. 
Thus, we restrict (u, v) so that 2 > Re (u + v) > |. 
Then, in the last expression for J(u,v;g) one may 
shift the j/-contour to (5 + e), encountering the pole 
at u + v — I with the residue 



R(u- 



C(2(2 



«)) 



<d(u,v;u + v~ §;#) 
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as well as those of the factor R(v + i)/£(l — 2v): 
we may assume, without loss of any generality, that 
u, v are such that all the residues in question are 
finite. This yields a meromorphic continuation of 
the continuous spectrum part, so that one may move 
(u, v) close to (^, ^) as far as Re (u + v) > 1 is sat- 
isfied; this condition is needed to have the last 6 
factor defined well. Then, shift the i^-contour back 
to the original. All the residual contribution com- 
ing from R(v+ ^)/C(l — 2i/) cancel out those arising 
from the previous shift of the contour. Only the 
pole at \-u-v contributes newly. The result- 
ing integral is regular at (5, 5); we get the factor 

Q A {u;g) = Q{\,\\v,g), 

9. Collecting all the above, we obtain 

Theorem. Under the above specifications, we have 
the spectral decomposition 

I(A; g) = M(A; 9) + 2Rc J ]T H v {\) Q A {V; g) 

I v 

where Ra = R, and M(A;g) is the value at (|,|) 
of the function 



R(u + v) 



■5(0) 



((2(u + v))» y ' C(2(2-u-«)) 
• + v — l)6(u, v; u + v — \;g) 

+ R(l - u — v)Q{u, v; § — u - v; g) 
+ R(u + v — l)6(u,u;u + w— §;g) 

+ R(l - u — v)Q(y, u: | — u - v; g) |. 



Albeit the Q factors in the last expression is defined 
so far only under the condition 2 > Re (u + v) > 1, 
the expression can in fact be continued to a neigh- 
bourhood of (i, |), for I(u,v;g) and all other parts 

in the spectral expansion of J(u,v;g) and J(v,u;g) 
are regular there. We could make the continuation 
procedure more explicit, using the property of T p 
mentioned in the first section, but the above suffices 
for our present aim. 

Our theorem should be compared with [3, The- 
orem] and [4, Theorem 4.2] which respectively deal 
with the mean square of automorphic i-functions 



associated with discrete series representations and 
with that of the product of two values of the Rie- 
mann zeta- function, i.e., the L- function associated 
with the continuous spectrum or Eisenstein series. 
Unlike those highly explicit results, admittedly the 
above asserts only the existence of a full spectral 
decomposition for I (A; g). The construction of the 
transform 6 a (•;<?) has to be made explicit in terms 
of the weight function g before one attempts any ap- 
plication; we shall take this task in our forthcoming 
works. Nevertheless, the present work could be a 
means to broaden the perspective of the theory of 
the mean values of the zeta and L-functions that is 
rendered in [7] via particular examples. 

Acknowledgement. We are greatly indebted 
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